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Abstract 



Despite the common belief, it is not always guaranteed that the light-front 
energy integration of the covariant Feynman amplitude automatically gener- 
ates the equivalent amplitude in light-front quantization. Our example of a 
light-front calculation with a fermion loop explicitly shows that the persistent 
end-point singularity in the nonvalence contribution to the bad component of 
the current, J - , leads to an infinitely different result from that obtained by the 
covariant Feynman calculation unless the divergence is properly subtracted. 
Ensuring the equivalence to the Feynman amplitude, we have identified the 
divergent term that needs to be removed from J~ . Only after this term is 
subtracted, the result is covariant and satisfies current conservation. The 
same calculation with the boson loop, however, doesn't exhibit such a sin- 
gular behavior and without any adjustment yields the result identical to the 
Feynman amplitude. Numerical estimates of the nonvalence contributions are 
presented both for the cases of fermion and boson constituents. 
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I. INTRODUCTION 



With the advent of Light-Front field theory, one can be quite hopeful to develop a con- 
nection between Quantum Chromodynamics (QCD) and the relativistic constituent quark 
model which is used in various electroweak form factor calculations. QCD provides a funda- 
mental description of hadronic and nuclear structure in terms of elementary quark and gluon 
degrees of freedom. It is very succesful in the perturbative regime, for instance when it is 
used in the explanation of the evolution of distribution functions in deep-inelastic scattering. 
There, the basic hard interaction is described by perturbative QCD and the Electro-weak 
interaction. The soft parts, distribution and fragmentation functions, are non-perturbative 
ingredients, which are less well understood. Lattice calculations are becoming increasingly 
accurate, but there is still much room and need for other non-perturbative approaches. 

It is part of the nature of the description of hadronic systems in terms of quarks and 
gluons that the characteristic momenta are of the same order or even very much larger than 
the masses of the particles involved. Therefore a relativistic treatment is called for. A very 
promising technique is Light-Front Dynamics (LFD), which treats relativistic many-body 
effects in a consistent way [|IJ. In LFD a Fock-space expansion of bound states is made. The 
wave function i/) n (xi, k^~, A«) describes the component with n constituents, with longitudinal 
momentum fraction Xj, perpendicular momentum and helicity Xi, i = 1, . . . ,n. It is 
the aim of LFD to determine those wave functions and use them in conjunction with hard 
scattering amplitudes to describe the properties of hadrons and their response to electroweak 
probes. 

Recently, important steps were taken towards a realization of this goal. In the work 
of BrodskyHiller and McCartor [[J, it is demonstrated how to solve the problem of renor- 
malizing light-front Hamiltonian theories while maintaining Lorentz symmetry and other 
symmetries. (The genesis of the work presented in may be found in || and additional 
examples including the use of LFD methods to solve the bound-state problems in field theory 
can be found in the review |T[). These results are indicative of the great potential of LFD 
for a fundamental description of non-perturbative effects in QCD. However, at present there 
are no realistic results available for wave functions of hadrons based on QCD alone. In order 
to calculate the response of hadrons to external probes, one might resort to the use of model 
wave functions. This way to estimate matrix elements was used by Ji et al. [Q. The same 
reasons that make LFD so attractive to solve bound-state problems in field theory make it 
also useful for a relativistic description of nuclear systems. Presently, it is realized that a 
parametrization of nuclear reactions in terms of non-relativistic wave functions must fail. 
LF methods have the advantage that they are formally similar to time-ordered many-body 
theories, yet provide relativistically invariant observables. 

Until now we have sketched a rather rosy picture for the application of LFD to hadron 
physics. However, not all is well and this is just the reason for the present investigation. 
Since the 1980's it has been assumed that the observables computed in the framework of LFD 
using the methods of perturbation theory are invariants, just as in covariant perturbation 
theory. Many authors have shown that LFD has this feature in particular cases and some 
years ago some general statements to the same effect could be made []5]|| . A case in point is 
the calculation of a current matrix element in quantum field theory. A typical amplitude is 
given by the triangle diagram. One encounters this diagram e.g. when computing the pion 
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form factor (see Fig. [I]). 



FIG. 1. Covariant triangle diagram 

The vertices denoted by \I/ are coupling constants in covariant perturbation theory. The 
hard scattering process is the absorption of a photon of momentum q by a (anti-) quark. 
In the LFD approach the covariant amplitude is replaced by a series of LF time-ordered 
diagrams. In the case of the triangle diagram they are depicted in Fig. ^. The first of these 



(a) 



(b) 




FIG. 2. Time-ordered triangle diagrams 



two diagrams is easily interpreted in terms of the LF wave functions \F However, the other 
diagram has a vertex that can again be written in the same way as before, but it contains 
also another vertex, denoted by that cannot be written as a LF wave function. It is a new 
element in LFD. We call this vertex the non-wave- function vertex. In order to compute the 
form factors in the time-like region, the contributions from these vertices must be included. 
Semileptonic meson decay processes also require the contributions from these vertices. One 
may try to avoid using them by choosing special kinematic conditions. It is known however 
that this will not be a simple task [|7|,|| . 

In the present work, we investigate the contributions from the non-wave-function ver- 
tices. We construct both the wave-function and non-wave-function vertices using pointlike 
covariant ones. The model used here is essentially an extension of Mankiewicz and Saw- 



icki's (1 + l)-dimensional quantum field theory model |IOfl , which was later reinvestigated 
by several others [§,|Il|-[l4| . The starting model wave function is the solution of the covari- 
ant Bethe-Salpeter equation in the ladder approximation with a relativistic version of the 

The covariant model wave function is a product of two free single 
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particle propagators, the overall momentum- conserving Dirac delta function, and a constant 
vertex function. Consequently, all our form factor calculations are various ways of evaluating 
the Feynman triangle diagram in quantum field theory. 

The importance of the contributions of the non-wave-function vertices was investigated 
in two cases: the electromagnetic form factors of a scalar and a pseudoscalar meson with 
spin-1/2 constituents. We also calculated the same form factor of a scalar meson with spin- 
constituents to see the spin effects. In 3+1 dimensions both the covariant and the LF 
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calculations are divergent and the model without any smeared vertex for the fermion loop 
is not well defined. This is in dramatic contrast to the case of spin-0 (boson) constituents, 
where regularization is not needed at all. In order to disentangle the issue of the non-wave- 
function vertices from the need of regularization, we performed our calculations in 1+1 
dimensions, where at least the covariant calculations for spin-1/2 constituents give finite 
results. 

It is commonly believed and widely used that the LF energy integration of the covari- 
ant Feynman amplitude generates the corresponding equivalent amplitude in the LFD. As 
we will show in this work, however, the equivalence between the LFD and the covariant 
Feynman calculation is not always guaranteed. The bad component of the current, J~, with 
spin-1/2 constituents exhibits a persistent end-point singularity in the contribution from the 
non- wave-function vertex. Unless the divergence in this contribution is properly subtracted, 
the singular behavior leads to an infinitely different result from that obtained by the co- 
variant Feynman calculation. Ensuring the equivalence to the Feynman amplitude, we have 
identified the divergent term that needs to be removed from J~ . Only after the identified 
term is subtracted, the result is covariant and satisfies the current conservation. 

In the next Section (Section |I|), we present both the covariant Feynman calculations and 
the LF calculations using the LF energy integration for the electromagnetic form factors of a 
pseudoscalar and a scalar meson with spin-1/2 constituents as well as the same form factor 
of a scalar meson with spin-0 constituents. Section |T| contains the numerical estimates of 
both the wave-function and non-wave-function vertices to the electromagnetic form factors 
of each case presented in Section O. The conclusion and discussion follow in Section |TV. 



II. CALCULATIONS 

The electromagnetic form factors can be extracted from the matrix elements of the 
current J M 

{p'\J»\p)=ie m {p'»+p»)F{q 2 ), (2.1) 

where e m is the charge of the meson and q 2 = (p' — p) 2 is the square of the four momentum 
transfer. If one uses the plus-component, J + = (J° + J 3 )/ \/2, the LF calculation gives two 
finite contributions, the wave-function part and the non- wave-function part, that add up 
to the covariant result, as expected. The importance of the non- wave-function contribution 
varies strongly with the momentum transfer and depends sensitively on the binding energy 
of the meson. For small values of q 2 and small binding energy, the wave-function part is 
dominant, but elsewhere the non- wave-function is essential for agreement between the LF 
calculation and the covariant results. 

The form factor can also be extracted from the minus- component of the current, J~ = 
( J° — J 3 ) / y/2. Covariance guarantees that it makes no difference whether the form factor is 
determined using the plus or the minus current matrix element. As LFD is not manifestly 
covariant, it may happen that J~ leads to a form factor different from the one determined 
using J + . As we show in this Section, the matrix element of J~ diverges in LFD. Unless 
one regulates J~ , the current cannot be conserved. To assure the current conservation, it is 
crucial to identify the term that causes the divergence. We have identified this term exactly 
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and found that it is an infinite function of the momentum transfer. If this infinite term is 
subtracted, the two LF contributions become finite as it must be in the conserved current. 
Moreover, their sum equals again the covaraint result as expected. However, the regularized 
LF contributions are different from the two parts of the form factor extracted from the plus 
current. The differences grow with increasing binding energy. 

A. Pseudoscalar Meson with the Fermion Loop 

The covariant fermion triangle-loop (Fig. [I]) for the pseudoscalar meson leads to the 
amplitude given by 

( '| J M | ) AN [ ^ ~ + p ■ p')^ + (k 2 — m 2 — k ■ p^p* 1 + (k 2 — m 2 — k ■ p)p >fl 
' J (2tt) 2 (k 2 — m? + ie)((k — p) 2 — m? + ie)((k — p') 2 — m? + ie) 

(2.2) 

where m is the fermion mass and N modulo the obvious charge factor e m is the normalization 
constant fixed by the unity of the form factor at zero momentum transfer. Even though we 
will present the unequal constituent mass case such as the kaon in the next Section of our 
numerical analysis, for the clarity of presentation we will focus in this Section on the equal 
mass case, such as the pion, only. 

The usual Feynman parametrization and the covariant integration yields 

,„,„,. , „ , IIS N f 1 , f 1 - x , 2m 2 (l - x) - (m 2 + (x + y - 1) 2 M 2 - xyq 2 )x , 

(PV» = z(p»+p'n- / dx / dy \ \ \— f- >- 2.3 

tt Jo Jo {(x + y)[x + y — 1)M Z + m z — xyq 1 ) 1 



where M is the meson mass. For q 2 = 0, the integration leads to fix the normalization as 

(2.4) 



1/N 1 '" 



7rM(4m 2 - M 2 ) 



Ml M 

j £Lrcti3iIi I 

4m 2 V4m 2 - M 2 WAm 2 - M 2 



In LFD, the form factor F(q 2 ) can be obtained by calculating either (p'\ J + \p) or (p'\J~\p). 
In principle, the result must be identical to the above covariant Feynman result regardless 
of which component of the current is used. However, this is not necessarily the case as we 
demonstrate in the following. 

First, the calculation of (p'\ J + \p) integrating out the LF energy k~ in Eq. (|2.2|) yields 
F(q 2 ) given by 

m N 



vr(2 + a) 



1 , (1 + a) 2 m 2 

dx-, 



o [m 2 -x(l -x)M 2 ][(l + a) 2 m 2 - (1 - x)(a + x)M 2 } 
(1 + a) (a - x)[x(a - x)M 2 - (1 + a) 2 m 2 } 
o """a[{a - x)(l + x)M 2 - (1 + a) 2 m 2 ][x(a - x)M 2 + (1 + a)m 2 } 

a 2 M 2 



+ dx 



(2.5) 

where a is given by q 2 = — f^f ■ In Eq. ( |2.5|) the first and second terms correspond to 
the contributions from the wave-function and non-wave-function vertices depicted in the 
first and second diagrams in Fig. |2], respectively. We have verified that the contribution 
from the non-wave-function part vanishes at a = 0, i.e. at q 2 = 0, indicating the absence 
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of a zero-mode contribution [D^] in the good component of the current J + . Adding both 
contributions in Eq. ( |2.5|) , we obtain 



2N{l + a)m 2 



naM[(2 + a) 2 m 2 - (1 + a)M 2 
2 + 2a + a 2 



4(l + a)m 2 + a 2 M 2 



=Artanh 



+ 




aM 



4(1 + a)m 2 + a 2 M 2 
M 



y/Am 2 - M 2 



: arctan 



V4m 2 - M 2 



(2.6) 



This result is identical to the form factor obtained by the covariant Feynman calculation 
given by Eq. ( |2.3|) and F(0)=1 gives the same normalization given in Eq. ( |2.4|) . 
On the other hand, the calculation of {p'\ J~\p) yields F(q 2 ) given by 



F{q 2 



N 



tt(2 + a 
(1 + a) 2 m 2 



1 + a)M 2 
dx 



dx 



'l-x) 



aM 2 



o [m 2 - x(l - x)M 2 }[(l + a) 2 m 2 - (1 - x)(a + x)M 2 } 
(1 + afm 2 + {(1 + a) 2 - (1 + x)}(a - x)M 2 



o a - x [(a - x)(l + x)M 2 - (1 + a) 2 m 2 ] [x(a - x)M 2 + (1 + a)m 2 ] 



where again it is apparent that the first and second terms correspond to the contributions 
from the wave-function and non-wave-function vertices, respectively. However, the non- 
wave-function part shows the end-point singularity coming from ^3^. Without subtracting 
the end-point singularity, the result is infinitely different from that obtained in the (p'\ J + \p) 
calculation. This is an astonishing result that deviates from the common belief in the 
equivalence of the LFD and the covariant Feynman calculation. Neither covariance nor 
current conservation is satisfied without a certain adjustment. In order to identify the term 
that must be subtracted, we rewrite the above equation as follows: 



N 



tt(2 + a) 



'l + a)M 2 [ dx- 
Jo 



'l-x) 



+ 



dx 



m 2 - x(l - x)M 2 }[(l + a) 2 m 2 - (1 - x)(a + x)M 2 } 
R(x, a 



a: 



x 



(2.8) 



where R(x, a) is defined by 

(1 + a) 2 m 2 [(l + afm 2 + {(1 + a) 2 



R(x, a) 



[I + x)}{a - x)M 2 



aM 2 [(a - x)(l + x)M 2 - (1 + a) 2 m 2 }[x{a - x)M 2 + (1 + a)m 2 } ' 



(2.9) 



In order to obtain the identical result to Eq. ( |2.6| ) from the (p'\J + \p) calculation, we find that 
R(a, a) must be subtracted from the numerator of the non- wave-function part integrand, 
i.e., 



F(q 2 



N 



n(2 + a) 



[1 + a)M z / dx 



[m 2 - x(l - x)M 2 ][(l + a) 2 m 2 - (1 - x)(a + x)M 2 ] 
R(x, a) — R(a, a) 



dx- 



a — x 



(2.10) 
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The subtracted term R(a, a) = ^| depends on the momentum transfer and never vanishes. 
While the subtracted result Eq. (|2.10f) is identical to Eq. (|2.6|) , it is interesting to note 
that the zero-mode contribution does not vanish in the (p'\J~\p) calculation as the non- 
wave-function part in Eq. fl2.10|) still survives even at q 2 = 0. However, the subtracted 
result Eq. ( |2.10| ) with the zero-mode contribution assures covariance and satisfies current 
conservation. 



B. Scalar Meson with the Fermion Loop 



The Feynman parametrization and the covariant integration of the fermion triangle-loop 
for the scalar meson gives the following amplitude: 



(p'\J^\p) = AN 



d 2 k (3m 2 + k 2 — p ■ p')k^ — (k 2 + m 2 — k ■ p')p^ — [k 2 + m 2 — k ■ p)p'^ 



(2k) 



(k 2 — m 2 + ie)((k — p) 2 — m 2 + ie)((k — p') 2 — m 2 + it) 



m 



„..N /"I f 1 -* , x((l-x-y) 2 M 2 
KF FJ nJo Jo y ((x + y){x + y -l)M 2 + m 2 - 

where the normalization N is again fixed by F(0) = 1 and given by 



xyq 



xyq 



2\2 ' 



(2.11) 



1/N 



Am 2 
7rM 3 



M 1 

~ ^Am 2 - M 2 



arctan 



M 



V 'Am 2 



M 2 



(2.12) 



In LFD, the calculation of (p'\J + \p) leads to F(q 2 ) given by 



F(q 2 



N 



dx-, 



+ 



(1 + a){2(l - x){2x + a) - (1 + a)}m 2 
tt(2 + a) [Jo [m 2 - x(l - x)M 2 ][(l + a) 2 m 2 - (1 - x)(a + x)M 2 } 
(1 + a) (a - x)[(l + a)(l + Ax - a)m 2 - x(a - x)M 2 ] 



dx 



o a[(a - x)(l + x)M 2 - (1 + a) 2 m 2 }[x{a - x)M 2 + (1 + a)m 2 } 



(2.13) 



where the contribution from the non-wave-function part again vanishes at q 2 = 0, indicating 



the absence of zero-mode contribution [13] in the good component of the current, J + . Adding 
both contributions, we find 



F(q 2 



2N(1 + a)m 2 



vraM 3 [(2 + a) 2 m 2 - (1 + a)M 2 ] 
8(1 + a)m 2 - (2 + 2a - a 2 )M 2 
4(1 + a)m 2 + a 2 M 2 



Artanh 
2a 



aM 



4(1 + a)m 2 + a 2 M 2 / 
M 



VAm 2 - M 2 



arctan 



VAm 2 - M 2 



(2.14) 



This result is identical to the form factor obtained by the covariant Feynman calculation 
given by Eq. ( [2.11 ) and F(0)=1 gives the same normalization presented in Eq. ( |2.11| ). 

However, the calculation of (p'\ J~\p) generates an end-point point singularity similar to 
the one observed in the pseudoscalar case. Defining the function 
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S(x, a) 



1 + a) 2 m 2 [(l + a)(l + 4x- 3a)m 2 + {a(2 - a) + (2a - l)x}(a - x)M 2 



aM 2 [(a - x)(l + x)M 2 - (1 + a) 2 m 2 ][x(a - x)M 2 + (1 + a)m 2 ] 



(2.15) 



we find F(q 2 ) given by 



F(g 2 ) 



N 



vr(2 + a) 



1 + a) fi, 4(1 + a)m 4 - 2(2 + a)(l - x)m 2 M 2 + (1 - x) 2 M 4 
ax 



M 2 Jo [m 2 -x{\- x)M 2 } [(1 + a) 2 m 2 - (1 - x) (a + x)M 2 ] 
S(x, a 



+ 



dx- 



a — x 



2.16) 



where the non- wave-function part again shows the end-point singularity coming from In 
order to obtain the identical result to Eq. (|2.14j) from the {p'\J + \p) calculation, we find that 
S(a, a) must be subtracted from the numerator of the non-wave-function part integrand, 
i.e., 



N 



ix(2 + a) 



'I + a) fi 4(1 + a)m 4 - 2(2 + a) (I - x)m 2 M 2 + (1 - x) 2 M A 
ax 



M 2 Jo [m 2 - x(l - x)M 2 ] [(1 + a) 2 m 2 - (1 - x) (a + x)M 2 } 
S(x, a) — S(a, a) 



+ / dx- 



a — x 



(2.17) 



The subtracted term S(a,a) = again cannot vanish. The zero-mode contribution is 

also visible since the second term in Eq. (|2.17| ) doesn't vanish even if q 2 = 0. The subtracted 
result Eq. ( |2.17| ), however, assures the covariance and satisfies the current conservation as 
in the previous calculation of the pseudoscalar meson. 



C. Scalar Meson with the Boson Loop 



For a comparison with the scalar constituents neglecting spin effects, we present in this 
subsection the calculation of F(q 2 ) for a scalar meson with a boson loop. The Feynman 
parametrization and the covariant integration of the boson triangle-loop for the scalar meson 
gives the following amplitude: 



(p>\J»\p) = N 



d 2 k 



(27r) 2 (k 2 — m 2 + ie)((k — p) 2 — m 2 + it)({k — p') 2 — m 2 + ie) 

2x - 1 



iV r 1 r l ~ x 

iip^ + p 11 ) — / dx dy— — 

KF F UnJo Jo y ( y (x + y)( y x + y-l)M 2 + m 2 



xyq 



2\2 ' 



(2.18) 



where the normalization N is given by 

1 



1/N 



2irM 2 (4m 2 - M 2 ) 



2(2m 2 — M 2 ) ( M 

—1-| = arctan . 

My/ Am 2 - M 2 \ ^m 2 - M 2 



[2.19) 



In LFD, the calculation of (p'\J + \p) leads to F(q 2 ) given by 



8 



N 



4?r(2 + a) 
+ / dx 



— 1 / dx 



[1 + a)(2x + a)(l - x) 



o [m? - x(l - x)M 2 ][(l + a) 2 m 2 - (1 - x)(a + x)M 2 } 
(1 + a) 2 (a - 2x)(a - x) 



o a[(a - x)(l + x)M 2 - (1 + a) 2 m 2 ][x(a - x)M 2 + (1 + a)m 2 ] 



^2.20) 



where the contribution from the non-wave-function part again vanishes at q 2 = 0, indicating 
the absence of a zero-mode contribution [13] in the good component of the current J + as in 
the fermion loop cases. Adding both contributions, we find 



F{q 2 



N(l + a) 



27raM 3 [(2 + a) 2 m 2 - (1 + a)M 2 



-J 4,(1 + a)m 2 + a 2 M 2 Artanh 



aM 



2a(2m 2 - M 2 ) 

H =^- arctan 

V4m 2 - M 2 



4(l + a)m 2 + a 2 M 2 , 
M 



V4 



M 2 



f2.21^ 



This result is identical to the form factor obtained by the covariant Feynman calculation 
given by Eq. ( |2.18 ) and F(0)=1 gives the same normalization presented in Eq. ( |2.19| ). 
Similarly, the calculation of (p'\J~\p) generates F(q 2 ) given by 



F{q 2 



N 



4vr(2 + a) 



a 



[l + a) 
M 2 

dx 



dx- 



2(1 + a)m 2 - (2 + a)(l - x)M 2 



m 2 - x(l - x)M 2 ] [(1 + a) 2 m 2 - (1 - x){a + x)M 2 ] 
2(1 + a)m 2 + a(a - x)M 2 



aM 2 Jo [(a-x)(l+x)M 2 -(l + a) 2 m 2 ][x(a-x)M 2 + (l + a)m 2 } 



. (2.22) 



Unlike the cases of fermion constituents, however, the non-wave-function part here does 
not exhibit the end-point singularity. Without any adjustment, we find that Eq. fl2.22p is 
identical to Eq. (|2.21 ) obtained from (p'\ J + \p). Thus, the result is automatically covariant 
and satisfies current conservation. However, the zero-mode contribution is still present in J~ 
current as one can easily see that the second term in Eq. ( 2.22|) doesn't vanish when a goes to 
zero. In the next section, we numerically estimate the importance of the non-wave-function 
part in all three cases that we presented in this section. 



III. NUMERICAL RESULTS 

We have estimated both the wave-function and non-wave-function vertices to the elec- 
tromagnetic form factors of each case presented in Section ||. For the three cases (A,B,C) 
presented in Section [TI|, we show the numerical results of the form factor calculated via the 
minus component as well as the plus component of the current. We denote the contributions 
from the wave-function and non- wave-function vertices as ^ and F^~\ respectively 
when the plus(minus) component of the current is used. The sum of the two contributions 
is denoted by F^~\i.e. F^ = F v + H + F+K 

For the numerical computation, we take for the experimental meson masses of the pion 
and the kaon m n = 0.140 GeV and rrix = 0.494 GeV resp. and vary the quark masses to in- 
vestigate the binding-energy dependence of the meson form factors. We call the pseudoscalar 
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meson with the equal quark masses and mass = 0.140 GeV the "pion". Likewise, the 
pseudoscalar meson with the unequal quark masses and the meson mass = 0.494 GeV 
is called the "kaon". 

The "pion" form factor with the quark mass m q = 0.250 GeV is shown in Fig. [3[ When 
the plus current is used, the valence contribution diminishes very quickly as Q 2 gets 
larger even though the normalization at Q 2 = is entirely coming from the valence part as we 
pointed out in Section |TJ. The crossing between and F^ appears at Q 2 below 0.05 GeV 2 . 
When the minus current is used, however, the valence contribution F^ is negligible even at 
Q 2 = and the entire result is essentially given by F~ v . The value of F^(0) corresponds to 
the zero-mode contribution in the minus current J~ and it is interesting to note that more 
than 90% of the form factor at Q 2 = is contributed by the zero-mode. Since F^ t = F t ~ t 
exactly coincide with the covariant result obtained by Eq.(3) for all Q 2 as they must, only a 
single solid line is depicted in Fig. The same applies to all of the other figures presented 
in this work. 

In Fig. |], we present the results for the "pion" by changing the quark masses in the 
following way: m q = 0.140, 0.077 and 0.0707 GeV, respectively. The closer m q is to m T /2 = 
0.07 GeV, the smaller the binding energy gets, and the slope of F^{~ at Q 2 = (or the 
charge radius) increases with decreasing quark mass, as expected. We find that the crossing 
between F^ and F^ occurs at a larger value of Q 2 and F^ becomes larger near Q 2 = 
as the binding gets weaker. This may be explained by the reduction of the probability to 
generate the non-wave-function vertex (or the higher Fock state) compared to the valence 
state as the interaction between the constituents gets weaker. Thus, in the weaker binding, 
F^ dominates over F^. Similarly, F^ becomes the main contribution near Q 2 = which 
is the only region where the form factor exists in the weak binding limit. Consequently, the 
zero-mode -F^(O) gets substantially diminished as shown in Fig. |j. 

The "kaon" form factor with m q = 0.25 GeV and m s = 0.37 GeV is shown in Fig. |5|. 
Compared to the "pion" case, the dominance of F^ extends to the larger Q 2 region and the 
crossover between Fj~ al and F^ is postponed beyond the range of Q 2 = 1 GeV 2 . The zero- 
mode i^(0) is also much smaller than in the "pion" case even though F~ v rises very quickly 
as Q 2 gets away from the zero range. As one can see in Figs. 5 (a) and (b), the contribution 
from the heavier quark struck by the photon is larger than that from the lighter quark struck 
by the photon. We have indeed confirmed that as m s gets larger only the contribution from 
the heavy quark struck by the photon dominates as expected. 

In Fig. |5], the form factors of the scalar partner to the "pion", which we call "s-pion" in 
the following, with m q = 0.25, 0.14, 0.077, and 0.0707 GeV are presented for comparison 
with the "pion" case. The basic features of F^ and F^ near Q 2 = are same as in the 
"pion" case because F^ al (0) = 1 must hold for any meson. However, as the binding gets 
weaker, we find that the "s-pion" form factors F^( change sign at smaller Q 2 -values. This 
indicates that electron scattering off the "s-pion" not only has zero cross section at a certain 
electron energy, but also that the electron energy that yields zero scattering gets smaller as 
the binding of the "s-pion" is weaker. Another dramatic difference from the pseudoscalar 
meson is the astonishing cancellation between F~ al and F~ v . Especially in the strong binding 
case, both F^ and F^ are huge but they cancel in a very remarkable way to yield exactly 
the same result as F^ t . 

In Fig. [F| the form factor of the scalar partner to the "kaon", i.e. "s-kaon", is plotted. 
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0.5 1 

Q ! 

FIG. 3. Pion form factor in LF calculation in 1+1 Dim. Pseudoscalar meson with spinor 
consitutents. M = 0.140 GeV, m q = 0.250 GeV. Fat lines correspond to the plus-current, thin 
lines to the minus-current. The solid line is the full form factor. It is the sum of the valence and 
the non-valence contributions. The separate contributions differ but the sums coincide. The form 
factor determined from the covariant amplitude is identical with the full form factor determined 
in the LF calculation. 



(a) (b) (c) 




FIG. 4. "Pion" form factor in LF calculation in 1+1 Dim. Pseusoscalar meson with spinor 
consitutents. (a) = 0.140, m q = 0.140. (b) m n = 0.140, m q = 0.077. (c) m n = 0.140, 
m q = 0.0707. The lines have the same meaning as in Fig. ||. Note the change in scale in panels (b) 
and (c). 
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FIG. 5. Kaon form factor. LFD calculation in 1+1 Dim. tuk = 0.494, m q = 0.250, m s = 0.370. 
(a) Charge +1 on the light quark, (b) charge +1 on the strange quark, (c) e q = 2/3, e s = 1/2. The 
lines have the same meaning as in Fig. |3[ 
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FIG. 6. "Pion" form factor in LF calculation in 1+1 Dim. Scalar meson with spinor con- 
stituents, (a) = 0.140, m g = 0.250. (b) m n = 0.140, m g = 0.140. (c) = 0.140, m q = 0.077. 
(d) m n = 0.140, m q = 0.0707. Note the change in scales in the latter two panels. The lines have 
the same meaning as in Fig. |3[ 
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FIG. 7. "S-Kaon" form factor LFD calculation in 1+1 Dim. Spinor quarks, mx = 0.494, 
m q = 0.250, m s = 0,370. (a) Charge +1 on the light quark, (b) charge +1 on the heavy quark, 
(c) e q = 2/3, e s = 1/2. The lines have the same meaning as in Fig. |3[ 



The basic feature is similar to the "s-pion". In Fig. |[ we show the results for the "s- 
pion" when the spinor quark is replaced by a bosonic quark. As we extensively discussed in 
Section [H], the subtraction of the end-point singularity is not required in the scalar quark 
case in contrast to the spinor quark case. In the scalar quark case, it is interesting to note 
that F~ al and reveal a large difference compare to the spinor quark case, while and 
are very similar to the spinor quark case. We find that the huge cancellation between F~ al 
and F~ v observed in the spinor quark case does occur in the scalar quark case only for very 
strong binding, and the most of F t ~ t is saturated by F^. However, the tiny contribution 
from F~ al near Q 2 = grows as the binding gets weaker and we have demonstrated the 
dominance of the valence part in the small binding limit regardless of the spin content, as 
discussed above. Fig. 9 shows the corresponding results for the "s-kaon" when the spinor 
quark is replaced by a bosonic quark. While the general features are similar to the spinor 
quark case, the bosonic quarks are more tightly bound together than the spinor quarks so 
that the charge radius of the meson is smaller than the case of spinor quarks as one might 
expect from the Pauli's exclusion principle. 
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FIG. 8. "Pion" form factor in LF calculation in 1+1 Dim. Scalar meson with boson con- 
stituents, (a) m n = 0.140, m q = 0.250. (b) m n = 0.140, m q = 0.140. (c) m w = 0.140, m q = 0.077. 
(d) m n = 0.140, m q = 0.0707. Note the change in scales in the latter two panels. The lines have 
the same meaning as in Fig. |3[ 
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FIG. 9. "S-kaon" form factor LFD calculation in 1+1 Dim. Boson quarks mj- = 0.494, 
m q = 0.250, m s = 0,370. (a) Charge +1 on the light quark, (b) charge +1 on the heavy quark, 
(c) e q = 2/3, e s = 1/2. The lines have the same meaning as in Fig. |3[ 



IV. CONCLUSION AND DISCUSSION 

In this paper, we have analyzed both the plus and minus components of the current 
quantized on the light-front to compute the electromagnetic form factors of pseudoscalar 
and scalar mesons. We considered spin- 1/2 consituents as well as spin-0 constituents and 
found dramatic differences between the two cases. Comparing with the covariant Feynman 
calculations, we notice that the common belief of equivalence between the manifestly covari- 
ant calculation and the LF calculation linked by the LF energy integration of the Feynman 
amplitude is not always realized. The minus component of the LF current generated by 
the fermion loop has a persistent end-point singularity that must be removed to assure co- 
variance and current conservation. A similar singularity was observed in the calculation of 
the fermion self-energy in || and |J. The plus component of the LF current, however, is 
immune to this disorder and provides a form factor identical to the one obtained doing the 
covariant Feynman calculation. This phenomenon is also associated with the spin-effect of 
the constituents because the calculation with the scalar(spin-O) constituents does not have 
the same symptom. Decomposing the LF amplitude into the wave-function and non-wave- 
function parts, it is interesting to note that the end-point singularity exists only in the 
non- wave-function vertex contribution. 
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Even after the singularity is removed, the minus component of the current sustains the 
zero-mode contribution while the plus component is free from the zero-mode. We have 
numerically estimated the importance of the non-wave-function vertices in all three cases 
that we discussed in Section |I[ We considered also the unequal constituent mass cases such 
as the kaon form factor. We find that the behaviors of F~ al and F~ v are trimendously different 
between pseudoscalar and scalar meson cases, while F^ al and F+ v have very similar features 
in both cases. The huge but remarkably exact cancellation between F~ al and F~ v shown in 
the scalar meson case persists even if the spinor quark is replaced by the bosonic quark. 
In the bosonic quark case, however, the binding between the constituents is stronger than 
the spinor quark case. We also notice that the zero- mode F~ v (0) diminishes as the binding 
gets weaker. Our results are quite consistent to the earlier observation fl3| exhibiting the 
smaller zero-mode contribution in the heavier quark systems. In all of these cases, our results 
show that if the meson is weakly bound then the contributions from the wave-function and 
the non-wave-function vertices to the plus current are separately almost the same as those 
for the minus current. Of course, their sums add up to the same number as the covariant 
Feynman result in both the plus and minus cases. 

The calculations carried out so far are semi-realistic as the model was 1+1- dimensional 
and only a point- vertex was considered. It is clear from a formal analysis of the 3+1- 
dimensional case, however, that a singularity of the same form will occur in the matrix 
element of J~ calculated in LFD regardless of dimensionality. A recent analysis of the 
Burkhardt-Cottingham sum-rule seems to reveal a similar divergence in the polarized spin- 



1/2 structure functions [ l5| . While the additional regularization may be provided by smear- 
ing the point- vertex with a realistic wave-function in the 3+1-dimensional covariant treat- 
ment of the current, the identification of the singular term as we achieved in this work would 
still be necessary for the smeared vertex cases. The importance of the non-wave-function 
parts may nevertheless differ numerically from the 1+1-dimensional case. This point is 
presently under investigation. 
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